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Abstract 

We study the two-boundary extension of a loop model — corresponding to the dense 
phase of the O(n) model, or to the Q — n 2 state Potts model — in the critical regime 
— 2 < n < 2. This model is defined on an annulus of aspect ratio r. Loops touching 
i-^h ■ the left, right, or both rims of the annulus are distinguished by arbitrary (real) weights 

which moreover depend on whether they wrap the periodic direction. Any value of 
these weights corresponds to a conformally invariant boundary condition. We obtain 
the exact seven-parameter partition function in the continuum limit, as a function of r, 
by a combination of algebraic and field theoretical arguments. As a specific application 
we derive some new crossing formulae for percolation clusters. 

^ ! 1 Introduction 

£Sj . The study of conformal boundary conditions (CBC) and boundary operators is one of 

the most fruitful aspects of the vast problem of solving two dimensional field theories 
and string theories. There are many reasons for this. In the equivalent 1+1 dimensional 
systems, CBC describe possible fixed points in quantum impurity problems, such as the 
multichannel Kondo problem [TJ, while boundary operators decide the stability of these 
fixed points as well as RG flows. In string theory, CBC describe possible branes, while 
RG flows in this language decide issues of (open string) tachyon decay f2j. In statistical 
mechanics, boundaries are roughly where couplings to the outside take place — for instance 
couplings to electrodes in quantum Hall effect type problems and their Chalker-Coddington 
type lattice formulations [3], [4] . 

From a more formal point of view, conformal field theories (CFTs) with boundaries are 
easier to tackle than their bulk counterparts when complicated features such as indecom- 
posability or non-unitarity are present. Most of the recent progress in our understanding of 
logarithmic CFTs for instance has come from the consideration of their boundary analogues 
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Taking a slightly different point of view, one of the basic objects in our understanding 
of CFTs has been the O(n) loop model, which led, in particular, to the development of 
deep links with the powerful SLE approach [8]. It is therefore no surprise that the issue 
of CBC for loop models should be a major problem. This issue has however been slow to 
evolve, in part for technical reasons: the Coulomb gas formalism, which is so successful in 
the bulk case, is very difficult to carry out in the presence of boundaries, for not entirely 
clear reasons [HE]. It took progress on the algebraic side — through the study of boundary 
algebras and spin models with general boundary fields — for the simplest families of CBC 
to even be identified properly. The works [TOj [11] finally showed that CBC were obtained 
in the dense loop model by simply giving to loops touching the boundary a fugacity m 
different from the one in the bulk. Associated conformal weights and spectra of conformal 
descendents were identified, and deep connections with the blob algebra [T2j |T3] (also called 
the One-Boundary Temperley-Lieb algebra) made. Subsequently, beautiful calculations in 
2D gravity [TU [15] recovered the results of [TOj [UJ. This will all be summarized in later 
sections. 

Our purpose in this paper is to continue the study of [TUJ [UJ and discuss situations with 
several boundaries and boundary conditions. In the case of calculations on an annulus for 
instance, this means giving different weights to loops touching the left, the right or both 
boundaries. We will end up in these cases with generating functions depending on seven 
parameters, and of course numerous potential applications to counting problems. 

Technically, the geometrical situation on the annulus has to do with understanding 
representations of Two-Boundary Temperley-Lieb algebras. We will devote a fair amount 
of time to this issue, which is essential in obtaining some of our results and conjectures. 
For early work and results in this direction see [TBI IT?] . 

The problem on the annulus is also deeply related with determining the spectra of XXZ 
hamiltonians with the most general boundary fields: this has been a very active question 
in the Bethe ansatz community lately [18] . We will in particular provide a complete answer 
for the spectrum of these hamiltonians in the scaling limit. 

More formally, the key question behind the calculations we will present is the determi- 
nation of fusion rules (and thus spectra of boundary conditions changing operators) in loop 
models. There are deep aspects to this, some of which will be discussed here but mostly in 
subsequent work. 

The paper is organized as follows. At the end of this introduction, we provide a summary 
of our results. Section 2 contains crucial algebraic preliminaries, where we define and study 
in particular the Two-Boundary Temperley-Lieb algebra. Section 3 contains Coulomb gas 
calculations where, thanks to a realization of the boundary algebras involving injection of 
charge on the boundaries, we are able to calculate a subset of all the critical exponents of 
interest. This is deeply related with the version of the problem involving XXZ chains with 
boundary fields that we also discuss briefly. Section 4 is the main section. Combining exact 
knowledge about hidden degeneracies (that come in part from the algebraic analysis in 
Section 2 — see also [33]), Coulomb gas arguments, and an educated guess on the structure 
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of boundary states, we are able to propose a formula for the most general, seven parameters 
dependent partition function. Section 5 contains various combinatorial applications, and a 
review of the few cases previsouly known, which our formulas all recover. In Section 6 we 
present a new combinatorial application, in the form of certain refined crossing formulae 
for critical percolation. Finally, Section 7 gives our conclusions 

Summary of the results: In this article we study a dense loop model on the annulus. 
Because of the boundaries and the non-trivial topology of the annulus, there are several 
types of loops, depending both on its homotopy (contractible or not) and which boundaries 
(none, only left, only right, or both) it touches. We distinguish all these kinds of loops by 
giving them different Boltzmann weights. For convenience we always ask the number of 
non-contractible lines to be even. This restriction will appear more clearly by defining the 
model on a lattice in the following section. 




Figure 1: The conformal loop model on the annulus. Different Boltzmann weights are given 
to the loops, depending on their topology (contractible or not) and if they touch a boundary. 
There can be a loop touching both boundariers if and only if there is no non-contractible 
loop (a). There is always an even number of non-contractible loops, and they are allowed 
to touch the boundaries (b). 

This model is endowed with conformal invariance, so we expect its partition function 
to be invariant under any conformal mapping. In particular we can study the model on 
a periodic strip of size L x N (L in the periodic direction), related to the annulus A = 
{z : Ri < \z\ < R 2 } by 

z ^ z = R 2 exp {i^j^j 
The geometry is characterized by the modular parameter 

q = (2) 
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where r = L/N = 2vr/log^. As a consequence of conformal invariance, the partition 
function must depend on the Boltzmann weights of the loops and on the modular parameter 
q only. It is a well-known result that the central charge of the dense loop gas is 



c = 1- 
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m(m + 1) 

where m is related to the Boltzmann weight of the bulk loops n through 

n = 2cos7, 7= -, m > 0. 

' ' ' m + 1 

Note that m is not restricted to be an integer. Let us also recall the Kac formula 



[(m + l)r — ms] — 1 
4m(m + 1) 



(3) 



(4) 



(5) 



Contractible 



Type 



Weight 



Parametrization 



Yes 
Yes 

Yes 

Yes 
No 
No 

No 



Bulk 

Boundary 1 

Boundary 2 

Both Boundaries 
Bulk 
Boundary 1 

Boundary 2 



n 
n\ 

ni2 
I 

h 



ni 



ri2 



n = 2 cos 7 
sin(n + 1)7 



sin H7 

sin(r 2 + 1)7 
sinr27 



ri G (0,m + 1) 



r 2 £ (0,m + 1) 



ni2 



sin(ri + r 2 + 1 - ri 2 )| sin(ri + r 2 + 1 + ri 2 )^ 
sinri7sin r 2 7 

/ = 2 cos x 

sin(ni + l)x 



h 



h = 



smuix 

sin(n 2 + l)x 
sin « 2 x 



Table 1: Loop weights and their parametrizations. 
Now we are ready to present the main result of this article. In full generality, the 
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partition function of the boundary loop model is given by 



-c/24 



z = TW^ q 

W n£Z 



q c / 24 v— v v— v sinfui + uo — 1 + 2j)x sin y h 



r 2 — 1 — 2n, — +r2 — 1+2 j 



P ( q "> j>in>o sinnixsinn 2 x 

g~ c / 24 sin(-ni + n 2 ~ 1 + 2j)xsinx h _ n+ 

P ^ j>\n>o sin-uixsinu 2 x 

I g~ C/24 y y sin(m - u 2 - 1 + 2j)xsinx /tri _ ro _ 1 _ 2 „. ri _ ro _ 1+2 . 
P(l) f^^ sinuixsin-n 2 x 

q-c/24 sin(— u\ — uo — 1 + 2j)Y sin y , N 

\ \ : t £ ! ±L± ±g' l -ri-r 2 -l-2n,-r 1 -r 2 -l + 2j fQ\ 

P ^ j>in>o sin-uixsin-^x 

where the seven parameters appearing are fixed by the seven different loop weights. The 
relations between all these parameters are given in Table [H Note that P(q) is our notation 

f°rrL>i(l-9 fc )- 



2 Some algebraic preliminaries 

Let us begin by introducing a few algebraic concepts that we will need throughout our 
discussion. Our model is the densely packed loop model on the tilted square lattice. A 
very convenient way to think about it is to view it as a face model (see Fig. [2]). Each face 
can be of two different kinds, corresponding to a horizontal or a vertical splitting of the 
loops. Each closed loop is given a Boltzmann weight n. The loops touching the boundaries 
are distinguished from the bulk ones in our model, and they are given different Boltzmann 
weights m, ri2 or n\2 if they touch the first boundary, the second one, or both of them. 
The total weight of a particular configuration is then n N n^ 1 n^ 2 n^ 2 where the TVj's are the 
numbers of loops of each kind. We shall later refine these weights to include information 
about the homotopy class (contractible or not) of each loop. 



2.1 The Temperley-Lieb algebra 

To begin with, we just drop the distinction of the boundary loops. Then partition func- 
tion of such a loop model can be reformulated in terms of local operators satisfying some 
commutation relations that will correctly count the closed loops. The trick is done by the 
celebrated Temperley-Lieb algebra [20], defined as follows. The Temperley-Lieb algebra 
TLn defined on N strands consists of all the words written with the N — 1 generators 
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Figure 2: A configuration of dense loops on the tilted square lattice. Loops touching at 
least once the first (resp. second) boundary are marked with a black (resp. white) blob. 



[1 < i < N — 1), subject to the relations 



(7a) 
(7b) 
(7c) 



The point of this definition originates in its graphic representation. Represent e, as an 
operator acting on N strands 

i i+1 



then (JZb)-((Zh) read respectively 



and 



i i+1 



i i+1 

KJ 

N 



N 



n 



i i+1 



i i+1 
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Each configuration of loops on a lattice of width N can be written as a particular word of 
the algebra TLpj (for example the configuration in Fig. [2l dropping the blobs coming from 
the boundaries, would be written eie^eTe^e2eie^e-je^e^e'i). In fact all the configurations 
can be generated by taking powers of the transfer matrix of the model, which reads 



(8) 





( N 






( N 


\ 


T' N = 


n 


1 + ei 




n 






[ »=i 

\i odd 


J 




\i even 


/ 



2.2 Boundary conditions and blob operators 

In the model we have just introduced, the loops touching the left or right boundaries of the 
lattice are not different from the other ones. We will refer to this particularly simple case 
as "free" boundary conditions. In this paper we deal with much more general boundary 
conditions. They consist in giving a different Boltzmann weight n\ (resp. n<i) to the loops 
which have touched at least once the boundary 1 (resp. 2). This is encoded in the transfer 
matrix by the addition of so-called "blob" operators b\ and 62 to the algebra TXjv- Their 
graphical representation consists of a black (resp. white) blob which marks the first (resp. 
last) strand. b\ acts as 



+ 



and 62 as 



They satisfy the defining relations 



N 



+ 



N 



i > 2 
bl 

eihei 



bi ei -- 
h 

niei 



(9a) 
(9b) 
(9c) 



and 



i < N - 2 
bi 

ejv-i&2ejv_i 



b2&i 
n 2 e N -i 



eib 2 



(10a) 
(10b) 
(10c) 



In what follows, we will assume that N is always even. In that case it is possible to have 
closed loops touching both boundaries. In order to count each of these loops with a weight 
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ni2, we impose the relation 

/ N \ I N \ I N 



i=i 



1 even 

which can be drawn as 




\ i=i / \ 

\i odd / \ 



KJ W 

J\f) 
'XfXJ 



i=l 
i even 



N 

"12 | ]j[ e, 

i=i 

\i even 



(11) 



ni2 n n • • • r\ r\- 



The generators ej, b\ and 62, subject to the relations (0, (f9l)- (fT0]l and the quotient (fTl]) . 
thus form the Two-Boundary Temperley-Lieb algebra on iV strands, denoted 2BTL^. A 
simpler case to which we shall sometimes refer is the One-Boundary Temperley-Lieb algebra 
IBTLn-, generated only by the e^'s and b\. The transfer matrix of the two-boundary loop 
model is then a generalization of Eq. j8|) 





( N \ 


f N \ 


T N = bib 2 


n i+e M 


n i+e A 




\/o~dd / 


\i even / 



It generates all the boundary loop configurations on a strip (see Fig. [2]) and gives the 
correct weights n to the closed loops in the bulk, and m, n-i or nu to the ones touching 
the boundaries. 



2.3 Generic irreducible representations of 2BTL 



Irreducible representations of the Temperley-Lieb algebra are well known, and are closely 
related to those of the quantum group SU (2) q . When q is not a root of unity 1 , the represen- 
tation theory of SU(2) q is essentially the same as the one of SU(2). In that case, the corre- 
sponding irreducible representations of the Temperley-Lieb algebra are said to be generic. 
The generic representations have a simple graphical interpretation, as the Temperley-Lieb 
algebra itself. The different modules (representation spaces) V s are given by configurations 
of (N — s)/2 half-loops and s strings. For example, consider the Temperley-Lieb algebra 
on 4 strands TL4, for which there are only three generic modules. 



V 



(WW) 



v 2 



v 4 



{MM} 



x Of course here q is not the modular parameter defined by (J2J. 
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The vertical lines are the strings, and the action of ej on two strings on the sites i, i + 1 is 
defined to be zero. 

Now if we work with the Two-Boundary Temperley-Lieb algebra 2BTLn (or with 
IBTLn), the generic representation theory is quite similar and has been studied in 
ITTT [T6] , The modules consist of the all states formed with half-loops and strings, but the 
half-loops can be marked with black or white blobs. Note that every black blob is neces- 
sarily on the left of every white blob. One can show [TTJ HU [19] that the dimension of Vo 
for IBTLn is 2^. For example, the module Vq is of dimension 16 for 2BTL^: 





V7 W , 


w v , 








y vj , 


w y , 



















This result will play an important role in the sequel, when we will have to deal with Coulomb 
gas arguments. 

Now consider the modules with strings. Half-loops between strings cannot be blobbed, 
since they are always separated from the boundary by at least one string, so b\ or b^ cannot 
act on them. The leftmost (resp. rightmost) string can carry a black (resp. white) blob. 
They can also be orthogonal to the blob, in the sense that they are not eigenstates of the 
projector b\, but of the orthogonal projector ("unblob") 1 — b\. Let us thus mark with a 
black (resp. white) square the action of 1 — b\ (resp. 1 — &2)- Then there is not only one 
module with s strings, but four, depending on the blob status (blobbed or unblobbed) of 
the leftmost and rightmost strings. For 2BTL± the modules with two strings are 

















H V> 








+ 4> V> 






• + W + 


» Vf" = < 


+ W t 


. vf = < 


+ W + 


► vr = < 


+ W + • 


W H 




VI + + 




W + + 






kJ + + 

V J 




kJ H 

V J 











Obviously these modules are related to each other by the blobbed/unblobbed transforma- 
tions 

61,2 -» 1 - 61,2 (13) 

and this symmetry between the projectors b\ and 62 and their orthogonals 1 — b\ and 1 — 62 
will indeed play some role in our analysis of the loop model. 
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2.4 Markov trace on 2BTL and the boundary loop model on the annulus 



Let us begin by dropping the blobs and the particular boundary weights, and start with 
the free/free partition function. Then the transfer matrix is given by JH]) in terms of the 
Temper ley- Lieb generators e^'s. There is an algebraic tool closely linked with the study 
of the Temper ley- Lieb algebra, namely the Markov trace, which is useful for our problem. 
We give only a naive definition of that tool here. The Markov trace of an element M of 
the Temperley-Lieb algebra is the number which is obtained by identifying the top and 
the bottom of the diagramatic representation of M, upon counting each closed loop with a 
weight n. For example, consider M = e±e2 on N = 4 strands: 



Although this object clearly depends on the number of strands N, we do not mention it 
explicitly. Given this naive definition, the Markov trace is associated with the Temperley- 
Lieb algebra itself, and does not require to know anything about its representations. Note 
that the Markov trace is not a trace in the common sense: it is not a sum over a basis of 
states of the diagonal action of M. However, there is a well-known relation between the 
Markov trace and the usual traces over the different generic modules of the Temperley-Lieb 
algebra. Using n = 2 cos 7, we have 



This relation calls for a few remarks. First, we use the convention that V s is empty if s > N 
or if s N mod 2. In particular, the sum is finite, and the terms depend on N. Then, note 



that sm s ( f n f l)7 is a polynomial in the variable n : = n, - n 2 _ ^ Eg±Z = n 3 - 2n, 



etc. They are actually Chebyshev polynomials of the second kind, U s (n/2). 

We consider now a modification of the Markov trace, which will be useful for our loop 
model. We can decide that we draw all the Temperley-Lieb diagrams in R — {0} instead of 
R and that when we compute the Markov trace, we give a weight n only to the contractible 
loops, and another weight I = 2cosx to the non-contractible ones. Again it is convenient 
to consider an example: 





(14) 




N N 
\ \ \ 
\ \ \ 



'III 
I I I 
/ / / 
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The modified Markov trace has quite the same structure as the previous one. In par- 
ticular, it can be decomposed on the usual traces over the generic modules exactly in the 
same way. 



This time the coefficient of each trace is a polynomial in I, and it is a remarkable fact that 
it does not depend on n at all. 

It turns out that we can define the Markov trace (or the modified Markov trace) in 
the same way for the Two-Boundary Temperley-Lieb algebra 2BTL, counting the blobbed 
loops with the appropriate weight rix, ri2 or n\2 (or l\, I2 for non-contractible loops in the 
case of the modified Markov trace). Recall that contractible loops touching both boundaries 
appear only if we work with an even number of strands N, so the number of strings must 
always be even, and we write s = 2j. Again, this object admits a decomposition on the 
usual traces over the different generic modules 



where the D^j are some polynomials in ni, ri2 and n (or l\, I2, I only if we are dealing with 

the modified Markov trace). The computation of the coefficients D%j can be achieved by 
various methods, see jTl] for a combinatorial proof or [19] for a more algebraic approach. 
The results are as follows. Let u\ and «2 be such that 




(15) 




a,(3=b,u 



(16) 



sin(ui + l)x 



(17) 



sinuix 



and 



sin(u 2 + Y)x 



(18) 



sin U2X 



then 




sm{u\ + U2 — 1 + 2j)xsinx 



(19a) 






sin-uixsin U2X 
sin(ni — U2 — 1 + 2j)x sin x 

sin U\X sin —U2X 
sin(— ni + u 2 - 1 + 2j)x sin x 

sin— mxsin U2X 
sin(— u\ — U2 — 1 + 2j)xsinx 



(19b) 



(19d) 



(19c) 



sm — mx sm —U2X 
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These equations are related by the blobbed/unblobbed transformation (JT3J) . To see this, 
note that the weight of a non-contractible loop marked with a black square (recall the black 
square stands for the action of 1 — 61) is simply I — l\ = & '° - n "m"^ 7 ■ The transformation 
(JT3J) has thus the effect of changing u\ into — m, or U2 into — u<i- These are indeed the 
transformations needed to pass from Ug to D%, or D%j, or DfJ. 

The physical interest of the Markov trace, or of the modified Markov trace, is that it 
counts automatically with the correct weight all the loops of a Temperley-Lieb element 
when the top and the bottom of the diagram are identified. This is exactly what we need 
to write down the partition function of our loop model. The transfer matrix on N strands 
is an element of the algebra 2BTLn, see (PT2l) , We want to work on an annulus of size 
L x N, so taking periodic boundary conditions in the L direction, the partition function of 
our loop model is just the modified Markov trace of a power of the transfer matrix. 

Z = Tr x Tk (20) 

Eq. (fT6j) gives the natural decomposition over the different modules 

Z = tr Vo Tk+ D 2jtr v ^Tk. (21) 

i>l 

a,(3=b,u 

This relation holds for every N and L. In particular, it must remain true in the limit 
L, N — > 00 with L/N fixed. Then if we introduce the (properly renormalized) characters 

K%? = { lim tr vaP T^\ (22) 



2j I L,N^oo V 2j 

\ J rcnorm. 

the conformal partition function will have the following structure 

Z = K + £ D$Kg. (23) 
i>i 

a,(3=b,u 

Hence, the computation of the conformal partition function has been reduced to the deter- 
mination of the characters K^f ■ 



3 Coulomb gas for the sector without strings 

In the previous section we explained why the partition function should have an algebraic 
structure coming from the Two-Boundary Temperley-Lieb algebra that we have just pre- 
sented, and hence can be decomposed on different sectors corresponding to the generic 
irreducible representations of 2BTL. Hence we are allowed to deal with each sector inde- 
pendently. This section is devoted to the computation by Coulomb gas arguments of the 
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character Kq, defined in the previous section as the trace over the module Vo of 2BTL^. 
First we detail how to obtain the parametrizations given in Table Q] in the Coulomb gas 
framework. These parametrizations have also a deeper algebraic origin associated with the 
Temper ley- Lieb algebra [HI HU US, US] , but the detailed discussion of this aspect will be 
deferred to H9l. 



3.1 A reminder: Coulomb gas on an infinite cylinder 

We now want to map our boundary loop model on a height model for which it is simpler 
to compute some quantities such as correlation functions or partition functions. We recall 
some classical arguments here. For the loop model on an infinite cylinder, the mapping is 
well-known. This must correspond to the limit N L (recall L is the periodic direction). 
First begin by giving each loop an orientation, then interpret the oriented loops as level 
lines for a height field h defined on the cylinder. The height varies by Ah = ±7r when 
upon crossing an oriented line. Each loop is counted with a weight e ±ry , depending on 
its orientation. The sum over the two orientations then gives the correct initial weight 
n = e i7 + e - * 7 to the original loop. Then it is generally argued that this model renormalizes 
to a free gaussian model with action 

S = -f- I {dhf d 2 x. (24) 



A a 



{(a + 7A) 2 -(7A) 2 }- (25) 



4vr _ 

This is however not sufficient to count correctly the loops which wrap around the cylinder. 
To do this, one has to add two charges e ±l ^ 1 ^ s>h at the ends of the cylinder. This modifies 
the scaling dimension of the vertex operator e tah to 

9 
4 

The value of g can then be fixed by the following argument. We started from a model in 
which the height difference when passing through a loop is Ah = ±7r, so the operator cos2/i 
should be marginal. This recquires A2 = 2 or A_2 = 2, so 

g = l±-. (26) 

7T 

The choice of the sign can actually lead to two different phases of the loop model, dense or 
dilute. We are working with a dense loop model, so we have to choose the solution g < 1. 
To finish, let us determine the central charge of this conformal field theory. The addition 
of charges at the ends of the cylinder has changed the behaviour of the partition function 
on the very long cylinder (JV S> L) by a factor e 7rN h/ 7r ) 7/9 . This is sufficient to identify the 
central charge, since we expect Z ~ e ~ ncN /® L [ n that limit, instead of e ~ nN / 6L without the 
addition of charges. Then we have 

(7/vr) 2 

9 

Defining m such that 7 = ^rj, this is nothing but the well-known formula jS]). 
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3.2 Boundaries in the height model 



Now we turn to the finite geometry of the annulus, and deal with the boundaries. Begin 
again by giving an orientation to each loop. The bulk loops are counted with a weight €~ %1 
if they are clockwise oriented, and e* 7 in the other case. The Temperley-Lieb generators 
ej's hence can be defined as acting on the orientations as shown in Fig. [3l It is not difficult 
to check that the e[s satisfy the relations ([7j also in the oriented loop language. 



Figure 3: Faces for the oriented loop model used for the Coulomb gas construction. The first 
line is just the identity in the Temperley-Lieb algebra, while the second line is a generator 
ej satisfying (J7J. 

We must find out how the blob operator b\ acts on the loop orientation. There are four 
different faces (triangles) with half oriented loops which can be combined to create b\ (see 
Fig. H]). Two of them conserve the orientation of the loop, which means that there is one 
arrow coming from the left side of the triangle, and one arrow entering it. The two others 
do not conserve it: both arrows point in the same direction. It is clear that the two faces 
which do not conserve the orientation cannot contribute to the weight of a loop touching 
only this boundary, because the orientation will be conserved everywhere else, so when we 
close the loop this contribution just vanishes. Now assume that the blob just adds some 
arbitrary phase factor e ±jri7 to a closed loop. Then requiring jib), the loop gets a weight 
m oc sin(r r + 1)7 instead of n = 2 cos 7. The correct normalization is fixed by jHfr). There 
remains one free parameter: the phase of the coefficients of the faces which do not conserve 
the orientation. We end up with the expression of the blob b\ given in Fig. [H where e* ri27 
is our free parameter. 

The same can be done for the second blob 62, so we actually have two free parameters 
coming from the boundary faces which do not conserve the orientation. Our problem has 
a global phase invariance, so one of them can be fixed, to give the expression of 62 shown 
in Fig. |H All the different loop weights can then be computed in terms of the parameters 
77, r2 and 772. The weights n\ and ri2 are given by the sum over both orientations 





sin(?7 + 1)7 



(28) 



m = 



sm?77 
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2i sin 777 



-17-17 



-Me 



-47*127 



-He' 



fl27 



1 

2i sin r27 



oir 27,' 



+ I 



-W27/ ■ + £ / 



Figure 4: Action of the blobs on the oriented loops. The orientation of the loops is not 
conserved by the blobs. 



and 



n 2 



sin(r 2 + 1)7 



(29) 



sin r27 

The weight of a loop touching both boundaries is a sum over four possible configurations 
of the orientations (see Fig. [5]), giving the parametrization 



sin 



fi + r 2 + 1 + r 12 



7 8111 



n + r 2 + 1 — r\2 



-1 



ni2 = — 

sin 777 sin r 2 7 

as claimed in the introduction (see Tabled]). 



(30) 



o 



-1 



4sinri7sin r 2 ^ 



e -t(ri+-r 2 +l)7^ ^ _g-iri27 
_|_ e i(n+r 2 +l)7 _ e iri27 



Figure 5: The four terms giving the parametrization (|30j) for the weight of a loop touching 
both boundaries. 



3.3 Spectrum in the sector without strings 

In the sector without strings, we can compute the conformal character Kq using Coulomb 
gas arguments. The previous prescription for the operators b\ and b 2 gives us a height model 
with the action (J24|) . with Neumann boundary conditions d y h(x,y = 0) = d y h(x,y = N) = 
0. Because of the boundary vertices which introduce some magnetic charge in the system, 
we see that there can be a difference of height if we turn once around the annulus: 

h(x + L,y) = 2pir + h(x,y), p G Z. 
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Clearly, p is the number of boundary vertices which inject charge in the system minus the 
number of those which take charge from it (see Fig. [4]). Such a configuration must be 
counted with a weight e* pri27 . 




x = L 



Figure 6: Coulomb gas on the annulus. The boundaries x = and x = L are identified, 
but there can be a difference of height h(x + L,y) = 2pir + h(x,y), p E Z because the 
charge is not conserved along a boundary (see also figure HJ. Here p = 3. 



In addition, we must treat the non-contractible loops. Note that a non-contractible 
loop which touch the boundary is no longer a loop in our prescription for the Coulomb gas, 
because it is broken in several half-loops on the boundary. However, the non-contractible 
loops which remain in the bulk must be counted correctly (see Fig. [7|. Remember that 
we do not want to compute the full partition function here, but only the character Kq 
corresponding to the representation of 2BTL without string Vo. Consider the following 
instructive example: we want to compute the trace over the module Vo of the following 
element of 2BTL 

eih = ^ | | 
Only 4 states in Vo contribute to this trace 



WW \JkJ WW wv>. 

It should be clear why there are exactly 4 = 2 2 states contributing to the trace. The top 
of the diagram corresponding to our element eib\ puts strong constraints on these states. 
More precisely, it gives all the information about the part which is disconnected from the 
bottom of the diagram [21]. Then if there are 2j lines (not blobbed, as in our example) 
going from the bottom to the top of the diagram, the states that contribute to the trace are 
exactly those we can form with half-loops, wathever their blob status is. In other words, 
the number of these states is the dimension of the module Vo on 2j strands, that is 2 2 K A 
more complete study of this is given in [TT|. This conclusion is sufficient for our discussion: 
each non-contractible loop in the bulk contributes with a weight 2 to the character Kq. 

This has an important consequence for our Coulomb gas construction : since each non- 
contractible loop that we cross when we go from one boundary to another is weighted by 
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2, we do not have to put some additional electric charge to correct their weight (unlike the 
case of the infinite cylinder that we discussed above). 



y = N 



y = 

Figure 7: The non-contractible loops in the Coulomb gas framework. Those touching the 
boundaries become an ensemble of half-loops between points on the boundary. Each non- 
contractible loop in the bulk contributes to the character Kq with a weight 2. 

Then, so far, we are able to count correctly the lines going from one boundary to another, 
and the non-contractible loops in the bulk. The next question is of course: what do we do 
with the vertices which conserve the charge, which should be given weights proportional 
to e"' 17 , e"" 27 , etc.? Our guess here is that they do not contribute to the universal part of 
the character, so Kq does not depend at all on r\ and T2- The reason for this is that the 
part involving r\ (resp. T2) in b\ (resp. 62) is diagonal, so it can be viewed equivalently as 
a field living on the boundary. We expect any such boundary field to flow towards a fixed 
boundary condition under RG, which should not depend on r\ (or T2). We have checked 
that conjecture numerically, by transfer-matrix diagonalization and extraction of the finite- 
size corrections. On the infinite strip of width N, the leading exponent h is related to 
finite-size corrections to the free energy per area unit through the well-known relation 

, _ , /boundary nh - TTc/24 f 1 \ 

~ /bulk + N + ^2 + °{WJ- (31) 

We have computed /at for sizes up to N = 18, then extracted the leading exponent h using 
(j31~T) up to order N~ 4 . Although we do not reach a very satisfying precision, our numerical 
results are compatible with the conjecture that h does only depend on r\2 (see Fig. [8]). 

Now we are ready to compute the character Kq itself. Let us decompose h(x, y) as 

h(x, y) = 2pn + h(x, y) 

where h{x + L,y) = h(x, y) and d y h(x, y = 0) = d y h(x, y = N) = 0. The integration over 
h gives the usual Z = g" 1/24 / 'P(q). Then we are left with the contribution of the height 
difference 2p-7r, counted with a weight e ipri27 as explained above 

Kq (X Zq 

pez P ez 
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Figure 8: Numerical results: we compute the largest eigenvalue of the transfer matrix, then 
exctract the leading exponent h from the finite-size corrections. Here we plot the quantity 
$ related to the exponent by h = 4^(^+1) versus r i ( a ) an d r i2 (b). Although the precision 
obtained here is not very satisfying, our conclusion is that $ (and hence h) does not depend 
at all on r\ and T2 (a), but we rather have = r±2 (b). 



where r = L/N. Now we want the expression of some Virasoro character, so we have 
to work with q = e~ nT , not e~ 2n / T . We perform the Fourier transform using the Poisson 
formula Yl p ~> Yin I dpe~ 2mnp . The sum becomes 

^ — ^ / ^ e -(i"9/-r)p 2 +ip(n27-27rn) 
n 

= f^l^ 1/2 ^ e -{^/Ag)(r 12 lh-2n) 2 
n 

= (T/g) l/2 ^2q h ^-^,r 12 -(c~i)/2^_ ^ 

n 

Normalizing the final expression such that the contribution of the identity operator (ri2 = 1) 
without its descendents is just q~ c ^ 24: , we end up with 

-c/24 

K ° = pTT q hri2 ~ 2n ' ri2 • (33) 

Note that, although this character depends only on r\2, it is not true that it does not depend 
on the loop weights ni and ri2, because all these parameters are linked by ([301) . Thus the 
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character Kq is a function of n, m, n-i and nyi as expected. Note also that, because of the 
parametrization (|30l) . we expect that Kq is invariant under 



ru -» -ru (34) 

and 

^12 — > H2 + 27r/7 = ri2 + 2(m + 1) (35) 
which is indeed the case for (|33j) . because of the symmetries of Kac's formula ((5j). 

3.4 Relation with the open XXZ spin chain 2 

The fact that the representation Vo of 2BTLjy is exactly of size 2 N and the oriented loop 
framework we developed above both suggest that there is some link with the celebrated 
spin 1/2 XXZ chain, with appropriate boundary conditions. We would like to develop a bit 
this subject in the following section. In fact, the equivalence between the representation Vo 
presented above and the so-called spin chain representation of 2BTLn was proved in [16J. 

It is well-known that the Temper ley- Lieb generators ej can be interpreted as a local 
Hamiltonian density, that is we can construct a simple Hamiltonian (here with the blob 
operators) 

N-l 

H = -Ai&i - X 2 b 2 - ei ( 36 ) 



i=l 



where Ai and A2 are two (so far unknown) constants. 



ei = -~ (oTof + o V A + cos 7 of of ) + (of - of +1 ) + ^ 

61 = — — — (sinsi7 of -I- COSS17 of + i cosri7 of) + ^ 

2smri7 2 

62 = 7—7- (sins 2 7 of + COSS27 of + i cosr 2 7 of) + - 

2 sin r27 2 



with 

ri2 = s 2 - si- (37) 

If we parametrize 

sin 7 sinri7 sin 7 sinr27 
Al = ~ ; : — 7 TT A2 = — ; ; — -, r-r. (38) 

sin <pi sm(ri7 + 0i) sin0 2 sm(r 2 7 + 4>2) 



2 This digression can be skipped at the first lecture. 
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then our Hamiltonian is, up to an irrelevant additive constant 



n 




i \ N ~ l 



(sin S17 erf + cos S17 of) + icotan^i cotan(ri7 + fa) <rf 



sin 0i sin(ri7 + fa) 

- — ; r4 TT ( sin s 27 <*n + cos s 2 7 cr^) + icotan^ 2 cotan(r 2 7 + fa) c>n 

sin02 sin(r 2 7 + <p2j 

Note that this is a Hamiltonian for the XXZ chain with non-diagonal boundary terms. 
This kind of Hamiltonian has been studied in great detail over the recent years [221 [T51 I23] . 
Note also that this Hamiltonian is not hermitian, which was already the case for closed 
boundaries with an SU(2) q symmetry [24] . 

Our derivation of the character Kq ([33]) applies directly to the Hamiltonian ([39]) . so we 
make the following conjecture about the spectrum of this spin chain. The universal part 
of the spectrum of TC does not depend on Ai and A 2 when these are positive real numbers. 
This has been discussed in some detail in [10] in the case of one boundary, and we expect 
this to be true also in the present case. Then the spectrum should depend neither on the 
parameters fa, fa, nor on 77 and r 2 . The only relevant parameter is the difference s 2 — s%, 
which is related to the weights of the loops touching both boundaries in the loop model via 
j37j and ([30]) . The spectrum of the XXZ Hamiltonian (|39j) is then given by 



N 

where vf = 7rs " 1 T is the "Fermi velocity" 



En = ~W ( h ^n,r 12 - c/24) . (40) 



4 The two-boundary partition function 
4.1 One-boundary case 



In [TO], it has been conjectured that the partition function on the annulus with one free 
boundary condition and one blob is 



q -c/24 



v ^sin(% + 2j)x JVi>ri+2j . ( ^sin(-ui + 2j) 



+ &ui i--"i^y _ ri ,_ ri+2j k (41) 



P (l) sinuix sin-ni7 

where we recognize again some polynomials in I = 2 cos x and l\ = h ™.^ 1 M ^ x . This partition 
function hence has the structure we have detailed in the previous section, but on the One- 
Boundary Temperley-Lieb algebra 1BTL. So far, we have failed to provide some Coulomb 
gas arguments to derive the exponents h rr . Strong numerical evidence has been given in 
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[TO] for general ri, and exact results have been obtained from Bethe ansatz when t\ is an 
integer [17]. Many of these results have been rederived since by Kostov using 2d quantum 
gravity techniques [TU [15]. Note the consistency with our computation of Kq from the 
previous section : the leading exponent we expect from (|33l) is h ri2jri2 . The one-boundary 
case should be recovered from ni = n and n\i = n\, that is T2 = 1 and ri2 = T\, We see 
that ^n,n is indeed the leading exponent appearing in ([4T1) . A more precise analysis of 
the relation between the character Kq for two boundaries and the character q hr i- r i /P(q) 
in ([4Tjl also exists, although it recquires more representation theory for the algebra 2BTL. 
We will report on this in [19J. 

4.2 Boundary states and the partition function 

Now we turn to the computation which is the core of this paper, and we determine com- 
pletely the partition function of our two-boundary loop model in the most general case. The 
main idea of this computation follows the work of Cardy on minimal theories [25]. We start 
from the one-boundary partition function Z\b, and compute its modular transform. The 
result is then interpreted as a scalar product between an initial boundary state \B\) and 
the final state (free), with an evolution operator ^o+Lo-c/12 = e -2nN/L(L +L -c/i2) } nsertec i 
(see Fig. [9]). Then we argue that this result together with the knowledge ofthe sector with- 
out strings is sufficient to guess the partition function of the form (i?2| q L o+L -c/i2 ^\). We 
conclude by computing the modular transform back, and get the general partition function 
in the form (|23p . 




Figure 9: Modular transform of the partition function. This corresponds to the "open" or 
"closed string channel" respectively. We have Zb 1 b 2 = q L o+ L o~c/i2 _ 
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Modular transform of the one-boundary partition function: We start from (|4T1) 
and use again the Poisson formula Ylj ~~ *■ SpJ* dje l27T ^ p , exactly as in (|32|) . 

Z 1B = (2g)- 1/2 rC/12 Y sin ("IX + ri(l/g)(p + xM) ~2a/4<,)[(x/«+p)*-hM*] (42) 
P(q 2 ) 4^ sinnix 

What with loops touching both boundaries? Something special must happen in the 
sector without strings, because of the loops touching both boundaries. The one-boundary 
partition function may be seen as a very special case of the two-boundary one, when n 2 = n 
and n\2 = n\. Thus, in the one-boundary partition function, the character Kq given by 
(f33]) is present, with the special value r 12 = r\. However, for a generic value of r 12 , the 
exponents h ri2 - 2 n,ri2 have no such special value. On the other hand, we expect all the 
exponents in the string sectors to be completely independent of n 2 , since they cannot 
depend on the weight of loops touching both boundaries. Hence, in this respect, the sector 
without strings decouples from all the other sectors. In particular, the formalism shown in 
Fig. [U should apply if we simply cancel the contribution Kq coming from the sector without 
strings. At the end of the computation, because of the form of the partition function (|23j) . 
it will be sufficient to add Kq with T\ 2 giving the correct weight n\ 2 to the loops touching 
both boundaries, see ([30]) . We have then 

ZlB - K (r 12 = n) = (2g)- 1/2 V ™ {UlX + ^/^(P^ + x)) ~2g/* a )[( x /*+p)*-hM 2 ] 

P(q 2 ) ^ sinmx 

-(^r 1/2 C^ E cos (^) <f {9/4) [ pM ^ )2 ] 

= (free\q Lo + Co - c l 12 \B{ Ul ,ri)). (43) 

Boundary states: Recall that the free boundary condition on the boundary 2 actually 
corresponds to u 2 = r 2 = 1. What we want to do now is to identify the terms of 



(B(u 2 ,r 2 )\q Lo+i °- c / 12 \B(u 1 ,r 1 )) = V (B(u 2 ,r 2 )\h a ) (h a \B( Ul , ri )) q 2h ° (44) 



where the sum runs over all the primary exponents appearing in (I43p . and the states \h a ) 
satisfy Lq \h a ) = Lq \h a ) = h a \h a ). We have to distinguish the two sets of exponents 
appearing in J43J). 



K = 1/4 S [(x/tt +p) 2 - (7 A) 2 ] : Eq. gHJ gives 

/rm iMi \ ii, r/ u (o ^-l/2 sin ("IX + ri(7/7rg)(p7T + x)) 
(-8(1,1) \h a ) {h a \B{ui,ri)) = (2g) 1 : . 

sinuix 
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so we can guess that the straightforward generalization holds 



(B(u 2 ,r 2 )\h a ) {h a \B(ui,n)) 
= (2gy 1/2 SlnX sin ( u iX + ri(j/Trg)(p7T +x))sm(u 2 x + r 2 (j/ng)(p7r 
sinnixsinn 2 x sin (x + (7/^) (f 71 " + x)) 

• = 5/4 [p 2 — (7/707) 2 ] : This time Eq. ([43]) seems to give simply 

(B(l, l)\h a ) (h a \B( Ul , ri )) = - (2g)- 1/2 cos (np 7 ) , 

a result which is independent of u\. This actually would lead to absurd conclusions. 
Indeed, we work in the string sectors, which all give non-contractible loops, so we 
expect all the terms to be affected somehow by the weights I, l\ and l 2 . The only 
terms which are completely independent of these weights appear in the sector without 
string, and we have cancelled this contribution. This contradiction comes from the 
fact that h a is even inp, so when we take the sum over all the exponents, only the even 
part of (B(l,l)\h a ) (h a \B(ui, ri)) remains. Inspired by the form of the coefficients 
we have already encountered, we can try the simple but non-trivial inclusion of the 
following odd term in p 

/rm 1 m s, \/u 10/ \\ (n \-V2 sin ( u ix) cos (npT) + cos(mx) sin (r lP j) 

{B(l, 1) \h a ) {h a \B(ui,ri)) = - (2g) 1 : , 

sin «ix 

leading to the generalization 

(B(u 2 ,r 2 )\h a ) (h a \B(ui,ri)) 

1/2 sinx sin (mx+prij) sin (u 2 x + pr 2 j) 
s\nu\xs\nu 2 x sm(x+P7) 

Modular transform: Although the (partly guessed) relations ([45]) and ([46]) seem quite 
complicated, they lead to quite a nice formula when we go back to the "open string channel" 
(see Fig. [9]). To see this, we need once again to perform a modular transform. The two 
sums appearing in (j4"4"l) are now 



(2g)- 1/2 ^ C/12 V sin x sin (u lX + pnj) sin {u 2X + pr 2 j) . 2(g/4) [ p 2_ (7/7rg)2 ] 
P(q 2 ) ^r^sinuixs'mu 2 x sin(x+P7) 



and 



Z _ = ( 2 g)-^i-^.y ^ 

P(q 2 ) sinnixsinn2X 

x sin (mx + n{-i/-Kg){p-K + x )) sin {u 2 x + r 2 ('j/Trg)(pTT + X )) ~2a/4 a )[(x/*+p) 2 -ft/nft 
sin (x + (l/irg) (pvr + x)) 
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We can compute the modular transform of each part independently, and add the contribu- 
tions in the end. Let us begin with Z + . The product can be decomposed as 

sin (u\x + prij) sin (u 2 + pr 2 j) 
sin (x + vi) 

= -Im ^ yZ eie2e i[(6iui+g2U2-l+2j)x+p7rg((£iri+e 2 r 2 -l)(7/7rg)-2j)] _ 
ei, 2 =±lj>l 

Z+ is then of the form 

z+= zZ z +^ e ^ 

ei,2=±l j>l 

and Z + (j, eip) is a sum over p 6 Z. Let us write i? = eiri + e2r2 — 1 and U = e\U\ + e 2 u 2 — 1. 
We use the Poisson formula Xlpez — * Ylnez / dpe t2npn to compute the modular transform 
of Z+(j,€i i2 ). Note that we also use <T 1/12 / 'P(q 2 ) = {t /2)- l / 2 q- l / 2A / P{q) as usual. The 
sum appearing in the computation is 

1 e ~(^/T)gp 2 e ipng(R-y/ng-2j) 



i2npn e -(Tr/T)gp 2 e ipng(R-y/ng-2j) 



g"-<si r l+ e 2 r 2 - 1 - 2n > e l'"l+ e 2 r 2 — l+2j _ 



£ 

= £ /* e< 

= ( T£ ,)l/2 q h R-'2n,R+2j-(c-l)/24 
n£Z 

Putting all things together, we get 

= 1 g~ c / 24 sin(eiMi + e 2 u 2 - 1 + 2j) x sin \ 3 

+ " 2 P( g ) ^ sme lUlX sme 2 u 2X 

(49) 

Now consider the case of Z_. The computation is quite similar. First we have to decompose 
the sinus product. Let x = (7/5) (p + x/ 71 "); then 

sin(rix + ui) sin(r2X + «2X) 
sin(z + x) 

= -Im £ E e i e 2 ei[(eiri+e2r2 ' 1_2n)x+(eiU1+e2U2_1 " 2n)xl 

ei i2 =±ln>0 

Then we use Poisson's formula X^pez — * Sjez / dpe t2lT ^ p to get in the end 

z _ = ±_<£^_ y y y + ^2 - 1 + 2j)xsinx ^ iri+MW .,. an , flfl+MW ,, +aj 

2 ei ^±iJ^^o sin ei«ix Sill € 2 «2X 
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and after some relabelling of the indices, we have 

1 q- c / 24 v - sin(eiui + e 2 u 2 - 1 + 2j)xsinx 

2 P(q) 4^ sinei-uixsine 2 n2X 



X 



g' lE l r l + E 2 r 2- 1 - 2n ' e l r l+ e 2 r 2- 1 + 2 J — ^ ^ g^ e l r l + ' ! 2 r 2- 1 ~ 2 "' e l r l+ e 2''2- 1 + 2 J \ . (50) 



n>0 n<0 



The two-boundary partition function: Adding the terms (|49l) and ([50]) . we find that 
the total contribution of all the string sectors is 

g-c/24 s i n ( eiUl _|_ e2U2 _ i _|_ 2j)xsinx , . , .,„. Ir .-s 

1 \ \ \ l_t t t t till llg' l e 1 r 1 +e2i-2-l-2n,e 1 r 1 +e2i'2-l + 2j 

If we now take into account the sector without strings and add its conformal character Kq, 
we obtain the partition function ([6]) of our loop model, as claimed in the introduction of 
this paper. This partition function has the form (|23l) as expected. We are now able to 
identify all the conformal characters corresponding to the different sectors. 

-c/24 

K h 2 ) = - ^ g ft ri+r2 -i-2n, ri+r2 -i+2 3 (- 52a ) 

-c/24 

n>0 

-c/24 

K%j = . ^2q h -ri+r 2 -l-2n,-r 1+ r 2 -l + 2j 

n>0 

-c/24 

n>0 

Note that these characters are all related by the blobbed/unblobbed transformation (fT3l) . 

5 Comparison with known results 

We would like to check our partition function against some known results from [25, 26, 9|. 
5.1 Critical percolation on the annulus 

As a first simple application of our result, we can turn to the critical percolation problem 
on an annulus. Critical percolation corresponds to I = l\ = I2 = n = n\ = 712 = n-12 = 1, 
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and in that case the partition function ((6)) is simply 



Z = l. (53) 

To see this, it is sufficient to note that our partition function reduces to the one-boundary 
partition function (|4"ll when l 2 = n 2 = n and n 12 = rii. The one-boundary case itself 
reduces to the free/free case when l\ = m = n. Then for I = n = 1, it is easy to see 
that Z = 1 using Euler's pentagonal identity Now if we want to compute, for example, the 
probability -P CT ossing that there is at least one contractible percolation cluster going from one 
boundary to the other, we have to vary the weight of loops touching both boundaries n\ 2 . 
Indeed, each percolation cluster is encircled by exactly one loop, and each cluster touches a 
boundary if and only if its surrounding loop touches it. Since we know that n\ 2 does only 
appear through r\ 2 in the conformal character Kq, we have 

Z(n 12 ) = l + K (r 12 )-K (r 12 = l) (54) 

and then 

^crossing = 1 — Z (n\ 2 = 0) 



K {r 12 = 1) - K (r l2 = 3) 
£ fc6 z ( q ^+k +q 6k^ +5 k+i_ 2q 



(55) 



which agrees with [26] 



5.2 Relation with Q-state Potts models 

It is a well-known result that the Q-state Potts model can be reformulated as a dense loop 
gas with the loop fugacity n = \J~Q. Let us recall here how this can be achieved. The 
Potts model can be defined on the square lattice, with a spin a x £ {1, . . . , Q} living on 
each site. Only the neighbouring sites interact, and the partition function is the sum over 
all the Potts spin configurations 

^Potts = J [ exp{K5(a x ,a x i)} (56) 

Potts <xx'> 

and this is rewritten as 

^otts =y, n ( i+ 5 ^ ,v ) ( 5? ) 

Potts <xx'> 

with v = e K — 1. Then the most important step is to interpret ([57|) as a random-cluster 
(Fortuin-Kasteleyn) partition function. For a given Potts configuration, the FK clusters 
live inside the Potts clusters. 

^Potts=E E ^ #FKb ° nds (58) 

Potts FKcPotts 
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Taking the trace over the Potts configurations gives the Potts partition function in its 
Fortuin-Kasteleyn representation 



Zp 



'otts 



E 

FK 



V 



#FK bondsQ#FK clusters 



(59) 



In the FK representation, the mapping to the loop model is obvious: one has to draw all 
the loops which encircle the FK clusters or the clusters on the dual lattice (see Fig. [TOj) . 
Now let iV be the number of loops, C the number of clusters and C* the number of dual 
clusters. Clearly, N = C + C* . Moreover, Euler's formula gives C = C* — #FK bonds + 
^lattice vertices. Then, up to an unimportant global factor, (|59l) becomes 



Zp, 



otts 



Loop 



E 4,' 



^ \ #FK bonds 



N 



(60) 



It is well-known that the Potts model is critical when it is satisfies the self-duality relation 
vj\[Q = 1. In that case, ([60]) is exactly the partition function of a loop gas with fugacity 

VQ- 




Figure 10: Mapping from the Potts model to the loop model. The black structures are the 
FK clusters, while dual clusters are in grey. 
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It is not difficult to generalize the previous discussion to the boundary case. Let us 
assume that the Potts spins living on the boundaries are restricted to some subsets S\ and 
S 2 C {1,...,Q}. Let 

Qi = \Si\ Q2 = \S 2 \ Qi2 = \S!nS 2 \ (61) 
then taking again the trace over Potts configurations we get the following relation instead 

of (USD 

Z Potts = £ v#™ h °^Q c Q?Q?Q& (62) 

FK 

where C is the number of bulk clusters, and Ci, C 2 , C12 are the number of FK clusters 
touching the boundary 1, 2, or both of them. Introducing the number of loops of the same 
type iV"i, N 2 , N12, it is clear that Ni = Ci, N 2 = C 2 and -/V 12 = C\ 2 since each boundary 
cluster is encircled by exactly one boundary loop. For the bulk loops, this is different because 
each one can encircle either a FK cluster or a dual cluster, so we still have N = C + C* . 
Now Euler's relation gives C + C\ + C 2 + C\ 2 = C* — #FK bonds + # lattice vertices. Up 
to a global factor, (|62l) becomes 

/ „ \ #™ bo„d, - , » „ , „„ 

^ v^i l^J (63) 

Again, we can impose the self-duality relation v = \[Q and then the identification of the 
loop weights is straightforward: 

At this point we have given the correct weight to all contractible loops. In this article 
we are interested in a loop model on an annulus, so we have to take care about the non- 
contractible loops. This turns out to be non-trivial, and rather crucial if we want to recover 
some known partition functions of the Potts model on the annulus. The subtlety comes from 
the non-contractible FK clusters which touch both boundaries, which must be restricted to 
the set Si H S 2 . However, in the loop model these configurations are those with exactly two 
non-contractible loop, each one touching one boundary. Such configurations are counted 
with a weight lil 2 7^ Qi 2 . To solve this problem, we must identify the term coming with 
the coefficient lil 2 in the loop partition function JB)), and give it the correct weight Qi 2 to 
get the Potts partition function. 

Let Zi x i 2 be this term in the loop partition function. To identify this term, it is necessary 
to analyse carefully the polynomials (fT9l) . These can actually be written in terms of the 
Chebyshev polynomials of the second kind U n (x), as pj] 



D b 2 ) = hl 2 U 2j - 2 {l/2) - (h + l 2 )U 2j ^{l/2) + t/ 2i _ 4 (Z/2) (65) 
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with similar expressions for the other polynomials D^j , obtained by using the blobbed/unblobbed 
transformation (fl~3l) . which maps l\ on I — l\ and/or I2 on I — ^2- With those relations the 
identification of is straightforward, noting that the constant coefficient of the polyno- 
mial U2n IS ( — 1)". 

^ = W - K% - K% + (66) 

Thus we have found the precise relation between our loop partition function and the Potts 
one 

Zpotts = -^loop + (Q12 — hh) %hh (67) 

where all the loop weights are given by jMJ for the contractible loops, and I = n, l\ = ni, 
?2 = n2 for non-contractible ones. Of course we could have improved slightly the mapping 
by distinguishing Potts clusters according to homotopy. However, this will not be necessary 
to recover the known results about the Potts model. 



5.2.1 Ising model 



We would like to use (|67|) to recover some results about the Ising model on an annulus, 
which appeared in [2H [9]. Assume for example that the Ising spins are fixed to + on the 
first boundary and to — on the second one. This corresponds in our formalism to Q = 2, 
Qi = Q2 = 1 and Q12 = 0. Then all the parameters (see Table Q] for the parametrizations) 
of the loop model are fixed: we have 7 = x = 71 /A, u\ = r\ = 112 = T2 = 2, ri2 = 3. Eq. j6|) 
then gives 



-c/24 



oop 



neZ 



P(q) ^ 2 



^3-2n,3+2j 



n>0 



-2 sin (-1 + 2j) ~ ^ q h ~ 1 ~ 2n -- 1+ ^ + sin (-5 + 2j) ~ ^ q h - 5 - 2n -- 5+2 ^ 1(68) 

n>0 n>0 J 

and adding the term (|66l) as in (|67l) . we get the Ising partition function 



Z 



+/- 



g -c/24 



neZ 



7 -c/24 



3-2n,3 



+ 



E 



P(q) £ 2 



V^sin(3 + 2j)~ + (-l) J ')2> 



r ?l3-2n,3+2j 



n>0 



-2 ( V2 sin (-1 + 2j) - + (-iyj ^ q h ~^-^ 

n>0 

+ (y/2 sin (-5 + 2j) | + (-1) J ) ^ q h-5-2n,- 5 



(69) 



n>0 



29 



Consider the second term between the brackets, which comes with a factor —2. Consider 
this term twice, and once make the reindexation j ^ j + 2, n — > n — 2, and the second time 
j ' ^ j — 2, n — > n + 2. The first term thus obtained cancels almost all the terms in the first 
sum between brackets, and the second one almost all those of the third sum. Collecting 
what remains after these cancellations, we have 

-c/24 

7 1 ^2 q h3 - 2n - 3 



1 „— c/24 - _ 



3>1 



2 P(g) 

^ (\/2sm (-5 + 2j) ^ + (-1)') (q h -^s+v + g fc -3,-s+2^ 



i>3 



2 (q h - 5 - 2n '- 3 + q h ~^A \ (70) 



n>0 

Now we write 2j = 8k + 2, 8A: + 4, ... and then 

-c/24 

7 , — Z \ " „^3-2n,3 



„-c/24 - 

_|_" J g fc 3,3+8(fc+l) _|_ qhl,3+&(k+l) _ q h 3,3+2+Sk _ ^l,3+2+8fc 

fc>0 

— g ft -l.-5+8(fc+l) — q h -3,-5+8(k+l) _|_ gh-l,-5+2+8(*+l) _|_ g^-3,-5+2+8(fc+l) | 

_-c/24 



P(q) 

w/ n>0 



sums 

-c/24 



= £ 1^^3,3+8*: _|_ q h l,3+Sk _ g^3,5+8fc _ g^l,5+8fe| ^2) 



^ ^-«-an,-3 + g ^-l-2n,3^ 
n>0 

Recall the Kac formula j5]) to see that h- r - s = h r>s , so all the terms combine to form the 

? -c/2< 

kez 

Again use the Kac formula and m = 3 (recalling that 7 = ^?pj) to see that /i3,3+8fc = ^3,5-8fc 
and then 

w fcez 

Here we recognize the Rocha-Caridi formula, and we conclude that 

Z+/- = Xi,3 (74) 
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as expected from [25]. We could have done the same calculation for the spins fixed to + on 

both boundaries. The computation is exactly as the previous one, except that Q12 = 1 so 
ri2 = 1 this time. This would have led to 



which is again a result of Cardy [25]. The other boundary conditions, such as free/+ for 
example, reduce to a computation with the one-boundary partition function, which has 
been studied in [10]. Again all the results agree with those of [25J. 

5.2.2 Three-states Potts model 

When Q = 3 the Potts spins have three colours A,B,C. For example, we can compute 
the partition function with all spins fixed to A or B with equal probability on the first 
boundary, and to B or C on the second one. We have then Qi = Q2 = 2, Q12 = 1, so the 
parameters (see Tabled]) of the loop model are 7 = 7r/6, u\ = r\ = 112 = = 2, r\2 = 3. 
The computation is exactly as in the Ising case. Eq. (|67l) gives 



z +/+ = Xi,l 



(75) 




c/24 



A3-2n,3+2j 



n>0 




(76) 
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Once again we see that the double sums actually collapse to give 



c/24 _ 

^3-27i,3 



Z AB/BC ~ -p^2^ q 



n£Z 



5^ (2sin(-5 + 2j)- + (-l)J') + q h -3,-wj 



,3+2j 



i>3 

-3 ^ ^-i-a-.a + q 

n>0 
,-c/24 



^ ^ ^g^l,3+12fc _ gh\ : -3+12k _|_ g^3,3+12fc _ g^3,-3+12fc \ {77) 



wy kez 

Using the Rocha-Caridi formula, we finally obtain 

Zab/bc = Xl,3 + X3,3 (78) 

which agrees with |25j . All the results from this reference concerning the Potts model can 
be deduced from our loop partition function j6]), with the relation (|67l) . 



6 Refined crossing formulae for percolation on the annulus 

It should be obvious that the seven-parameter partition function j6]) harbours many more 
geometrical applications than the known ones presented in the preceding section. As an 
illustration we present here just one simple example. 

Consider the continuum limit of critical percolation on an annulus of aspect ratio r = 
L/N, and recall that q = e _7rT . Let Pq, be the probability that no cluster wraps the periodic 
direction, and let P^ be the probability that there are precisely j > 1 wrapping clusters 
which are moreover constrained by the values of the indices a, (3. When a = b (resp. 
a = u) the leftmost cluster is constrained to touching (resp. to not touching) the left rim; 
f3 similarly constrains the behaviour of the rightmost cluster. 
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Since Z = 1 we have obviously 

Pa = z(x = |,»i 



a/3 _ 



a,(3 



pbb _ 



puu 
3 



(2j)! \d lX J 
1 

(2^2;! \i),\J 
1 



= l,«2 = l) 

= 1, u 2 = 1 



2j ~ 2 d Ul d U2 Z(ui = l,ix 2 = 1) 



(duiM (^2^) 



(2j)! V^X 



9 \ 2i 

' » Z(«i = -1,1X2 = -1) 



(79) 



and since Pj w = P^ b by symmetry, this suffices determine all P^ '. Note also that by an 
easy duality argument we have = P"" for j > 1. 

We find the following explicit results for j < 3, here given to order ~ g 8 : 

1 4 7 q 10 a 13 c 16 K 

P = 1 - gs - gs + 2g 2 - 2gs + 2^ - 2q~i + 4g 4 - 4gir + 4<f - 5g^ + 8<f 



if 



p«fe 

puu 



pub 



r>uu 
*3 



19 7 22 o 

-8gT + 8q 7 - 10g^ + 14g 8 H 

1 4 Q 7 Q 10 13 c 16 fi 

gs - 2g + gs - 2g 2 + 2gs - 4g 3 + 6g^~ - 6g 4 + 8g3" - 12g 5 + I3q^ - 16g 6 
+20qf - 28q 7 + 30gf - 38g 8 + •■■ 

13 K 16 R 19 v 22 

g-r + 4g 5 - 2gir + 2g b - 3gi?" + 6g 7 - 5g^~ + 7q H 



g-g 



gs -g 

10 



= 9 + q - 2g^ + 2g 4 - 2g^" + g 5 - 4g^~ + 3g 6 - 6gir + 10g 7 - 10g^ + 12g 8 + 

10 13 c 16 R 19 7 22 a 

= q~ +q~ - 2g 5 + 2g~ - 2g b + 3g~ - 7g 7 + 5g~ — 9g H 



g 5 + g 6 + g 8 + 
g 7 + g 8 + • • • 



(80) 



and P^ u = g 3 + • • • . The evaluation of the complete series for aspect ratio r = 1 leads to 
the following numerical values: 



j 


^a,/3 ^7 


pbb 

3 


pub pbu 

3 3 


puu 

3 





0.6364540018880 








1 


0.3615910259567 


0.2770671481561 


0.0413139498152 


0.0018959781702 


2 


0.0019548143402 


0.0018959781702 


0.0000293394720 


0.0000001572261 


3 


0.0000001578149 


0.0000001572261 


0.0000000002943 


0.0000000000002 



These values could presumably be verified by numerical simulations in a square geometry. 
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7 Conclusion 



In this article we have studied a densely packed loop model on the annulus, with general 
loop weights that distinguish the two boundaries and the homotopy class of the loops. The 
main result is the exact seven-parameter continuum limit partition function j6j). We have 
verified that a range of special cases of this expression agree with existing results in the 
literature, and used it to derive new refined crossing probabilities in critical percolation. 

The directions for future work are quite numerous [19]. Let us discuss briefly a few of 
them: 

• Distinguishing both rims of the annulus by non-trivial boundary conditions is related 
with properties of 1BTL boundary condition changing operators. Resulting fusion 
rules are encoded in the result J6|). An intriguing — and to our knowledge novel — 
feature is that the fusion here depends on a parameter n\2 which is unrelated to those 
characterizing the two individual 1BTL operators. 

• Specializing the two-boundary model to simpler cases gives rise to a rich hierarchy of 
restrictions. For instance, the two-boundary model with n\2 = n\ and Ui = n becomes 
the one-boundary model, and with n\ = n this in turn becomes the standard ("zero- 
boundary") Temperley-Lieb model. Moreover, in each case there are "magical" values 
of the weights, typically corresponding to one of the r-type parameters taking an 
integer value. Each of these restrictions corresponds to the disappearence of some of 
the states in the transfer matrix, the vanishing of certain eigenvalue amplitudes, and 
the reorganization of the Hilbert space into new modules. There is a rich algebraic 
meaning of this truncation hierarchy. 

• The present work pertains to the dense phase of the O(n) model. In the dilute case 
the possibilities are richer: in addition to the boundary-specific n-type weights, one 
can weigh differently the boundary monomers depending on the type of loop to which 
they belong. This gives rise to several surface transitions. Some of those will be 
insensitive to the values of the r-type parameters, others will correspond to the usual 
swapping of indices (i.e., h r>s — > h s ^ r in the Kac formula), and yet others lead to 
genuinely new behaviour. 
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